Roughening of close-packed singular surfaces 
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An upper bound to the roughening temperature of a close-packed singular surface, fee Al (111), is 
obtained via free energy calculations based on thermodynamic integration using the embedded-atom 
interaction model. Roughening of Al (111) is predicted to occur at around 890 K, well below bulk 
melting (933 K), and it should therefore be observable, save for possible kinetic hindering. 

PACS: 68.35.Bs, 68.35.Md, 68.35.Rh 



Roughening [Q-H is one of the most fundamental 
phase transitions at surfaces, yet probably the most elu- 
sive. The roughening of vicinal surfaces |^,^ is gener- 
ally accepted to be a transition of infinite order of the 
Kosterlitz-Thouless class. The extremely weak free- 
energy divergence at the critical point implies frustrat- 
ingly slow variations in space and time of whatever order 
parameter is chosen to characterize the transition. This 
makes predictions on roughening a challenge for atom- 
istic simulations techniques, this being not the last of the 
reasons why statistical mechanics models have tradi- 
tionally been the dominant approach to this problem. 

The roughening of singular faces poses additional prob- 
lems. Vicinal surfaces roughen as the (mostly configu- 
rational) entropic free energy related to step meander- 
ing prevails over the cost of step and kink formation; 
on vicinals, where steps already exist by construction, 
this occurs generally at temperatures well below melting. 
Singular-face roughening, on the other hand, requires 
step formation to begin with. Singular faces, therefore, 
roughen at much higher temperatures, so much so that 
roughening is thought to be preempted by melting in 
most cases, especially on close-packed faces. 

Here we use a simple approach to predict the roughen- 
ing transition temperature of a singular surface, based on 
free energy calculations performed with an atomic-level 
finite-temperature simulation technique (the embedded 
atom method coupled with Monte Carlo thermodynamic 
integration). We calculate the free energies of several 
vicinals to the singular face, and estimate the tempera- 
tures at which the free energy of each vicinal becomes 
lower than the singular. Since roughening is phenomeno- 
logically identified with the appearance of hills and val- 
leys of arbitrary height on the surface, we assume that 
roughening will be fully developed at the temperature at 
which the steepest and most costly vicinal will be favored 
over the low-index face. To obtain an internally consis- 
tent and low-error-bar estimate, we calculate the crossing 
temperatures of the free energies of several vicinals with 
progressively shorter terraces, with the free energy of the 
singular surface; we then obtain Tr as the extrapolated 
crossing temperature of the shortest/most costly vicinal. 



To be definite, here we estimate an upper bound to Tr 
for any Al surface, and find it to be ~ 890 K, well below 
the bulk melting temperature of 933 K. 

To obtain such upper bound, we study Al (111), which 
is expected to have the highest roughening temperature 
among the low-index faces, being the most closely packed. 
Also, it is stable ||] up to the bulk melting temperature, 
and predicted to sustain overheating 0. The vicinals of 
Al (111) we consider here are Al (8 8 10), Al (557) and 
Al (335), obtained by miscut of the (111) plane at an 
angle of ~ 1, 9, and 14 degrees respectively. There exist 
two kind of steps on Al (111), namely the Ill-facetted 
and the 100-facetted. The latter are energetically more 
costly, and our vicinals belong to this second class. In 
the notation of Lang [|l^ , bearing out directly the inter- 
step distance, these faces are denoted as [9(111) x (100)], 
[6(lll)x(100)] and [4(lll)x (100)], respectively, meaning 
(say) 6 rows of a (111) face separated by a (lOO)-faceted 
step. These vicinals lay on the (Ill)-(IOO) line of the 
stereographic map of the fee lattice The steepest 
vicinal on this line is Al (113), or [2(111) x (100)] : its 
appearance should set the occurrence of fully developed 
roughening. Here we first simulate straight-step vicinals, 
and then estimate the correction due to kink formation 
by simulating one kinked vicinal. 

Free energies are calculated via the embedded atom 
method and thermodynamic integration. The embedded 
atom method [^ is a fairly reliable method to predict 
structural and thermal properties of metals. Its main ad- 
vantage is its moderate computational cost, and ensuing 
high numerical accuracy achievable within the method's 
bounds. The disadvantages are essentially that the choice 
of materials to be simulated is restricted by the availabil- 
ity of accurate potentials (constructing which is a science 
in itself), and that the embedded atom method, by its na- 
ture of effective interatomic potential, is not as accurate 
as first principles methods. This inherent inaccuracy is 
attenuated for Al by the highly refined parameterization 
of Ercolessi and Adams [|l^ , built to reproduce a large 
database of ab initio energy and force calculations. Re- 
cently [|ll| the Ercolessi-Adams model has been further 
refined to cure minor inaccuracies in the description of 
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surface diffusion and high-energy scattering. 

Thermodynamic integration is adopted because the 
roughening transition occurs (if at all) well above the De- 
bye temperature 400 K for Al bulk), and it is therefore 
imperative to properly include anharmonic effects in the 
free energy of the relevant surfaces. While useful at lower 
temperatures, the commonly adopted quasi-harmonic ap- 
proximation is not very reliable at high temperature, as 
shown by recent simulations on Al (100). In ther- 
modynamic integration p^ , the potential energy of the 
system is progressively switched on, through a parameter 
A, starting from a reference system whose free energy is 
known: 



V(A) =AW-(l-A)Uh 



(1) 



with W and U^. the potentials of the actual system and 
of an harmonic crystal. Since fl^ , 
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the free energy at a fixed temperature Tj-cf is 

F„f = Fa=i = Fa=o + / (Uh - W)a dA , 
Jo 
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The integrand is calculated by Metropolis canonical 
Monte Carlo, and the A=G value is known by construc- 
tion. By the thermodynamical free energy-enthalpy rela- 
tion 
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the free energy in the interval [TrcfjT] is 



F(T) = T 
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The integrand is calculated again by canonical Monte 
Carlo. The surface free energy per unit area is 
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[Fsiab (T)-NFbuik(T)] 



(6) 



where Fbuik (T) is the bulk free energy per atom, Fgiab 
(T) is the free energy of the N-atom simulation slab, and 
A its surface area. Whenever appropriate, thermal ex- 
pansion is accounted for with an expansion coefficient 
given by the ratio a(T)/a(G) of the lattice constants at 
temperatures T and zero, obtained by an NPT simu- 
lations |11|. Free energies are calculated in supercells 
containing 450 to 600 atoms depending on the orienta- 
tion. Each Metropolis Monte Carlo run was ^3x10^ 
steps long. We estimate error bars of 0.5% in the sur- 
face free energy, and about ±10 K in the crossing points 
and Tr. 
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FIG. 1. Free energy temperature dependence for Al (111) 
and vicinals. 

In Fig. |l] we report the free energy vs temperature for 
Al (111) and vicinals. At first, Al (111) is favored. At 
higher temperatures, vicinals with progressively shorther 
terraces become favored free energy-wise. The crossing 
points are T=773 K for Al (8 8 10), T=875 K for Al (557), 
and T=914 K for Al (335). 

The crossing points in Fig.|l| tend to "accumulate" to- 
wards a finite value as the terraces become progressively 
shorter. This suggests to identify Tr with the "accumu- 
lation point" of this sequence. To quantify it, we fit a 
polynomial through the crossing points just obtained as 
function of interstep distance, and define Tr as the tem- 
perature value corresponding to the interstep distance on 
the vicinal surface with the shortest terrace within our 
class of (lOO)-faceted, straight-stepped vicinals, namely 
Al (113), whereby the interstep distance is ^ 5.5 A. At 
T^Tr as just defined, all the vicinals (within our re- 
stricted class) are favored over Al (111), so that arbi- 
trarily large and composite fluctuation can appear in the 
surface profile. The result is displayed by the upper curve 
in Fig. 1^: the roughening temperature estimate is 930 K, 
very close to the melting temperature of Al bulk (theo- 
retical §: 939 K; experimental [|: 933 K). 

To refine the prediction, we first note that only un- 
detected straight steps have been considered so far. On 
the other hand, at finite temperature kinks will form on 
steps. Kinks affect the free energy of the stepped sur- 
face both indirectly because their very existence allows 
step meandering, and directly via their formation internal 
energy, and vibrational entropy due to their vibrational 
modes. We consider that the latter free-energy variation 
will be captured accurately by a simulation. We neglect 
instead the step meandering-related configurational en- 
tropy, based on previous work on vicinal surfaces p^ ] 
suggesting that the configurational contribution is negli- 
gible compared to the vibrational below the roughening 
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transition. 
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FIG. 2. Crossing points of surface free energies of vici- 
nals and singular surface. Upper curve: straight steps; lower 
curve: kinked steps. Tr is defined as the temperature corre- 
sponding to the Al (113) interstep distance. 




FIG. 3. Top view of kinked Al (557) as studied in free en- 
ergy calculation. As all other cells, it contains two periodi- 
cally-repeated steps per side. 

The simulation of vicinal surfaces with kinks is de- 
manding in periodic boundary conditions; here we re- 
strict to a single case, kinked Al (557), chosen because 
of its favorable geometry. Each side of the simulation 
slab, depicted in Fig. |^, contains one straight and one 
kinked step. The latter exhibits two kinks, with a rel- 
atively low linear density of 0.05 A^^. The number of 
atoms is preserved by this procedure, as required by nu- 
merical considerations. As shown in Fig. ^, the kinked 
Al (557) turns out to have a crossing point with Al (111) 
at T=845 K, with a reduction of 4% over the straight- 
step value. Assuming that the other crossing points are 
lowered by about the same amount due to kinks, and ap- 



plying the same procedure as before, we find Tr= 887 K 
(lower-laying curve in Fig. ^). This is a strong upper 
bound because accounting for lower-cost (lll)-faceted 
steps should lower this figure. In addition, account for 
meandering will also lower (moderately) our estimate. 

Roughening has not been reported for any (111) face 
so far. The predicted Tr is rather close to, but lower 
than the melting temperature, so it is quite conceivable 
that roughening of Al (111) could be observed. Our pre- 
diction concerns energetics, however. Kinetic effects are 
not considered in any way. However, Al (111) was ob- 
served Q in Medium Energy Ion Scattering experiments 
to remain stable up to the melting temperature. Also, 
molecular dynamics simulations showed Al (111) to 
be stable for at least 2 ns up to 1088 K, or 150 K above 
bulk melting. While the length and time scales accessible 
in simulation are not comparable with those of relevance 
in roughening, this is an indication that kinetics may 
play a role, slowing down or hindering the transforma- 
tion. Thus, it is possible that experiments aiming at the 
observation of the roughening of Al (111) predicted here, 
may have to observe the surface over time spans of hours, 
or produce "nucleation" defects by e.g. nanoindentation. 
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FIG. 4. Lowering of the free energy crossing point of Al 
(557) with Al (111) due to the presence of kinks. 

As a further check on the predictions based on the 
embedded atom Al potential, reinforcing the plausibility 
of our estimate, we calculate Tr for vicinals within the 
Terrace-Ledge-Kink (TLK) model of Villain et al. 
through the relation 

K^ Wm ^Wo/kBTR._ 

keTR 

Here Wm is the energy needed to move a step by one row 
towards a neighboring step m-l-1 atomic rows away, and 
Wq is the kink formation energy. This expression is valid 
for W,n < T < Wo, which is the case here. The value 
of K depends upon the details of the underlying theory, 
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and it equals 2 for the original TLK model; values of 2 
for Cu (113) ^ and 2.1 for Ag (115) |] have been sug- 
gested based on experiments or MonteCarlo simulations 
on vicinals. We evaluate these parameters from total en- 
ergy calculations on slabs containing at least 5 steps per 
slab side, and comprising from 1700 to 4000 atoms de- 
pending on orientation. The parameter Wm is calculated 
removing one complete atomic row of step-edge atoms. If 
N is the total number of atoms, and L that of step-edge 
atoms, the total energy for row remowal is 



L W,„ — En-l — [En — LEb], 



(8) 



with En-l and En the internal energy of system after and 
before row removal, and Eb the bulk energy per atom. 
Wm is thus defined per atom. For a kink we only remove 
half a row, creating two kinks: 



2Wo = E^ 



[En - (L/2) Eb] - W„ 



(9) 



^N-L/2 

with En-l/2 the internal energy of the slab after half-row 
removal. 

For Al (335) we find Wm=3 meV, Wo=112 meV, Tr= 
411 K; for Al (557) we find W^^l meV, Wo=108 meV, 
Tr= 314 K; for Al (8 8 10) we find Wm=0.1 meV, 
Wo=106 meV, and Tr= 209 K. These values are quite 
comparable with results of previous investigations on 
stepped metal surfaces [ p7| , p^ . Our numbers for Al (335) 
are compatible with those inferred from STM measure- 
ments on Ag (115) [^, which has the same step-step 
separation: Wni= 3 meV, Wo=114 meV, and Tr=427 
K. [The (115) face consists of (lll)-faceted steps sepa- 
rated by a (100) terrace four atomic rows wide, whereas 
the (335, has (100) steps and (111) terraces.] Concerning 
Tr of Al (111), den Nijs et al. observed [^ roughening 
of Ni (115) at about 200 K, and estimated 420 K for the 
roughening of Ni (100), the nearest singular face on the 
stereographic plot. Our value of 412 K for Al (335) simi- 
larly suggests that our upper bound of 890 K for the asso- 
ciated singular (111) is quite plausible. Our predictions 
for both singular and vicinal faces await experimental 
verification. 

In summary we have calculated an upper bound to 
the roughening temperature of a singular metal surface 
using an atomistic simulation method. Our results for 
Al (111) suggest that roughening may occur apprecia- 
bly below melting, and be therefore observable, save for 
kineting hindering. 
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